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Abstract 

Social foraging behaviour of Escherichia coli bacteria has 
recently been explored to develop a novel algorithm for 
optimization and control. One of the major driving forces of 
Bacterial Foraging Optimization Algorithm (BFOA) is the 
chemotactic movement of a virtual bacterium that models a 
trial solution of the optimization problem. However, during 
the process of chemotaxis, the BFOA depends on random 
search directions which may lead to delay in reaching the 
global solution. This paper comes up with a hybrid approach 
involving Particle Swarm Optimization (PSO) and BFOA 
algorithm called Bacterial Swarm Optimization (BSO) for 
designing Static Synchronous Series Compensator (SSSC) in 
a power system. In BSO, the search directions of tumble 
behaviour for each bacterium are oriented by the 
individual's best location and the global best location of PSO. 
The proposed hybrid algorithm has been extensively 
compared with BFOA and PSO. Simulation results have 
shown the validity of the proposed BSO in tuning SSSC 
compared with BFOA and PSO. Moreover, the results are 
presented to demonstrate the effectiveness of the proposed 
controller to improve the power system stability over a wide 
range of loading conditions. 
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Introduction 

The power transfer in an integrated power system is 
constrained by transient stability, voltage stability and 
small signal stability. These constraints limit a full 
utilization of available transmission corridors. Flexible 
AC Transmission System (FACTS) is the technology 
that provides the needed corrections of the 
transmission functionality in order to fully utilize the 
existing transmission facilities and hence, minimizing 
the gap between the stability limit and thermal limit 
(Kundur, 1994). 

Recently, there has been a surge of interest in the 
development and use of FACTS controllers in power 
transmission systems (Lee, and Sun, 2002; Sen, 1998; 



Ngamroo, 2001; Ngamroo, and Kongprawechnon, 2003; 
Zhang, 2003). These controllers utilize power 
electronics devices to provide more flexibility to AC 
power systems. The theory and the modeling 
technique of SSSC device using an Electromagnetic 
Transients Program (EMTP) simulation package is 
presented in (Sen, 1998). SSSC is designed in 
(Ngamroo, 2001) to stabilize the frequency of 
oscillation in an interconnected power system. This 
SSSC located in series with the tie line between any 
interconnected areas, is applicable to stabilize the area 
frequency of oscillations by high speed control of tie 
line power through the interconnections. A robust 
design of the lead/lag controller equipped with the 
SSSC for stabilization of frequency oscillations is 
discussed in (Kongprawechnon, 2003). A multi control 
functional model of SSSC for power system analysis is 
described in (Zhang, 2003). This model can be used for 
steady state control of one of the following parameters: 
1) the active power flow on the transmission line; 2) 
the reactive power flow on the transmission line; 3) the 
voltage at the bus; 4) the impedance (precisely 
reactance) of the transmission line. A robust clamping 
controller based on Fuzzy Logic Controller (FLC) is 
introduced in (Chen, Lie, and Vilathgamuwa, 2004). 
The only input signal for this damping controller is the 
real power measurement at the location of the SSSC to 
generate the modulation index for controlling the 
injected voltage of the Voltage Source Converter (VSC) 
while its phase angle is required to remain constant 
with respect to local reference voltage vector. A new 
Genetic Algorithm (GA) based approach for optimal 
selection of the SSSC damping controller parameters in 
order to shift the closed loop eigenvalues toward the 
desired stability region is designed in (Kazemi, 
Ladjevardi, and Masoum, 2005). The dynamic 
operation of both Static Synchronous Compensator 
(STATCOM) and SSSC based on a new model 
comprising full 48-pulse Gate Turn Off (GTO) VSC is 
investigated in (El-Moursi, and Sharaf, 2006). These 
models combined reactive power compensation and 
voltage stabilization of the electric grid network. The 
rate of dissipation of transient energy is used as a 
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measure of system damping in (Haque, 2006). This 
concept is applied to determine the additional 
damping provided by a STATCOM and SSSC. GA 
optimization technique to design FACTS based 
damping controllers for Single Machine Infinite Bus 
(SMIB) is applied in (Panda, and Padhy, 2007). The 
analytical expressions for this additional damping are 
derived and compared for classical model of a simple 
power system. The influence of SSSC and STATCOM 
on the synchronizing power and damping power of a 
SMIB is introduced in (Jowder, 2007). The impacts of 
different SSSC control modes on a small signal and 
transient stability of a power system is discussed in 
(Castro, Ayres, Da-Costa, and Da-Silva, 2007). The 
application of a SSSC controller to improve the 
transient stability performance of a power system is 
thoroughly investigated in (Panda, and Padhy, 2007). 
The performance of different input signals to the 
Power Oscillation Damping (POD) controller is also 
assessed. The transient energy is used as a tool to 
assess the effectiveness of FACTS devices to damp 
power system oscillations in (Murali, and Rajaram, 
2009). A systematic procedure for modeling, 
simulation and optimally tuning the parameters of a 
SSSC controller for power system stability 
enhancement is presented in (Panda, 2010). 

Several optimization techniques have been adopted to 
solve a variety of engineering problems in the past 
decade. GA has attracted the attention in the field of 
controller parameter optimization. Although GA is 
very satisfactory in finding global or near global 
optimal result of the problem; it needs a very long run 
time that may be several minutes or even several hours 
depending on the size of the system under study. 
Moreover swarming strategies in bird flocking and fish 
schooling are used in the PSO and introduced in 
(Kennedy and Eberhart, 1995). However, PSO suffers 
from the partial optimism, which causes the less exact 
at the regulation of its speed and the direction. Also, 
the algorithm cannot work out the problems of 
scattering and optimization (Rini, Shamsuddin, and 
Yuhaniz, 2011; Selvi, and Umarani, 2010). In addition, 
the algorithm pains from slow convergence in refined 
search stage, weak local search ability and algorithm 
may lead to possible entrapment in local minimum 
solutions. A relatively newer evolutionary 
computation algorithm, called BF scheme has been 
addressed by (Passino, 2002; Mishra, 2005; Fogel, 1995) 
and further established recently by (Ali and Abd- 
Elazim, 2011). The BF algorithm depends on random 
search directions which may lead to delay in reaching 
the global solution. A new algorithm BF oriented by 
PSO is developed that combine the above mentioned 
optimization algorithms (Korani, 2008; Biswas, 
Dasgupta, Das, and Abraham, 2007; Abd-Elazim, and 
Ali, 2012). This combination aims to make use of PSO 
ability to exchange social information and BF ability in 



finding a new solution by elimination and dispersal. 
This new hybrid algorithm called Bacterial Swarm 
Optimization (BSO) is adopted in this paper to solve 
the above mentioned problems and drawbacks. 

This paper proposes a new optimization algorithm 
known as BSO for optimal designing of the SSSC to 
damp power system oscillations. The performance of 
BSO has been compared with these of PSO and BFOA 
in tuning the SSSC damping controller parameters. 
The design problem of the proposed controller is 
formulated as an optimization problem and BSO is 
employed to search for optimal controller parameters. 
By minimizing the time domain objective function, in 
which the deviations in the speed, DC voltage and 
transmission line power are involved; stability 
performance of the system is improved. Simulation 
results assure the effectiveness of the proposed 
controller in providing good damping characteristic to 
system oscillations over a wide range of loading 
conditions. Also, these results validate the superiority 
of the proposed method in tuning controller compared 
with BFOA and PSO. 

Power System Modelling 

SSSC is installed in series with transmission line as 
shown in Fig. 1. The generator is represented by the 
third order model that comprising of the 
electromechanical swing equations and the generator 
internal voltage equation. The IEEE type ST1 excitation 
system is used (Kundur, 1994). Details of system data 
are given in appendix. 

S = co B {co-\) (1) 

a = —\P -P -D(co-l)) ( 2 ) 
r v m e * 

j 

where P 1 and P are the input and output powers 

of the generator, respectively; tj and D are the inertia 

constant and damping coefficient, respectively; 8 and 
CO are the rotor angle and speed, respectively; ft>g is 

the synchronous speed. 




Ox: 
X: 

FIG. 1 SMIB WITH SSSC 
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The output power of the generator can be expressed in 
terms of the d axis and q axis components of the 
armature current and terminal voltage as following: 



P =Vji'j +v i 
e da q q 



(3) 



The internal voltage, E equation is shown below: 



E' = — -E' +—E~j + 
1 t', 1 t\ f d 



'do 



L do 



X d X d 
T 'do 



(4) 



where E ^ is the field voltage; is the open circuit 

field time constant; X ^ and X ^ are the d axis 

reactance and d axis transient reactance of the 
generator, respectively. 

Modeling ofSMIB with SSSC 

The SSSC is a VSC connected in series with the 
transmission line at its midpoint through an insertion 
transformer as shown in Fig. 1. The SSSC output 
voltage is defined by the following equation (Castro, 
Ayres, Da-Costa, and Da-Silva, 2007; Mandour, and 
Abd-Elazeem, 2009): 

= cVqq {cosy/ + j'siny/ - ) (5) 

Where c is the amplitude modulation ratio, y/ is the 
phase angle modulation ratio of the SSSC andV^^ is 
the SSSC DC voltage. 

The transmission line current is described by the 
following equation: 



I tL =I tLq + j I > 



tLd 



(6) 



The SSSC DC voltage differential equation is given 
below: 



DC c 



DC 



{I tLqCOW + I tLd sm ¥ } 



(7) 



The induced AC system voltage due to SSSC voltage is 
described by the following equation: 



bt 



■v h + jx h i 



b tL 



(8) 



Substitute from equations (5) and (6) into equation (8), 
and then divide it into d and q axis as following: 



V btq + J V btd = CV DC (C ° S V + J S ' m + 

The machine terminal voltage is described as following: 



V t = ^ X tL + X LB^L +V bt +V B= V q + jV d 



(10) 



From equations (6), and (9) the transmission line 
current in d and q axis are defined as following: 



tLq ^ X tL +X LB^ 



Wd- V btd 



- Vg sin S} 



'tLd 



7 q +V bt q +V B™V 



(11) 
(12) 



AC Voltage Regulator 

The AC voltage regulator controls (Bamasak, and 
Abido, 2005) the reactive power exchange with the 
power system as shown in Fig. 2. 

Where Kp , and Ki are the PI controller gains for 

t ac ac b 

AC voltage regulator, K , and 7^ to T, are the lead lag 

controller gains for additional controller in the AC 
voltage regulator circuit, U ^ is the output signal of 

additional controller in the AC voltage regulator 
circuit, and T^y , is the washout time constants for AC 

voltage regulator and its additional controller 
respectively. 
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FIG. 2 SSSC DYNAMIC MODEL OF AC VOLTAGE REGULATOR 

DC Voltage Regulator 

The DC voltage regulator controls (Bamasak, and 
Abido, 2005) the DC voltage across the DC capacitor of 
the SSSC controller as shown in Fig. 3. 
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FIG. 3 SSSC DYNAMIC MODEL OF DC VOLTAGE REGULATOR 

Where Kp , , and Ki^ c are the PI controller gains for 
DC voltage regulator, Vj^q is the SSSC DC voltage, 
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and r, is the washout time constant for DC voltage 

regulator respectively. To reduce the computational 
burden in this study, the value of the wash out time 
constants 7 ' and 7^, are fixed to 8 second, the values 

of and 7^ are kept constant at a reasonable value 

of0.05 second. The parameters of the AC, DC, and 
additional controller are to be determined via various 
optimization techniques. 

Objective Function 

In the present study, an integral time absolute error of 
the speed deviations, DC voltage and transmission line 
power of SSSC is taken as the objective function 
expressed as follows: 



siwi i ) 
J= \ t\\Aa)\ + \AV DC \+\AP t \)dt 





(13) 



The great advantage of this objective function is that 
the effect of SSSC signal is taken into consideration. 
The problem constraints are the SSSC controller 
parameter bounds. Therefore, the design problem can 
be formulated as the following optimization problem: 



Minimize J (13) Subject to 



kpT^p 



dc 



< Kp 



dc 



Ki™ m <Ki, <Ki™ ax 
dc dc dc 



Kp mm <Kp < Kp max 
' ac 1 ac 1 ac 



Ki mm < Ki < Ki max 
ac ac ac 



K mm < K < K max 



T mm <T <T max 



(14) 
(15) 
(16) 
(17) 



< r < t~ 



(19) 
(20) 



Hybrid BFOA-PSO Optimization Algorithm 

Overview ofBFOA 

Natural selection tends to eliminate animals with poor 
foraging strategies and favour the propagation of 
genes of those animals that have successful foraging 
strategies since they are more likely to enjoy 
reproductive success. After many generations, poor 



foraging strategies are either eliminated or shaped into 
good ones. The Escherichia coli bacteria that are 
present in human intestine also undergo a foraging 
strategy. The control system of these bacteria that 
dictates how foraging should proceed can be 
subdivided into four sections namely Chemotaxis, 
Swarming, Reproduction and Elimination and 
Dispersal (Mishra, 2005; Fogel, 1995). 

Chemotaxis 

The characteristics of movement of bacteria in search 
of food can be defined in two ways, i.e. swimming and 
tumbling together known as chemotaxis. A bacterium 
is said to be 'swimming' if it moves in a predefined 
direction, and 'tumbling' if moving in an altogether 
different direction. Mathematically, tumble of any 
bacterium can be represented by a unit length of 
random direction <p(j) multiplied by step length of that 
bacterium C(z'). In case of swimming, this random 
length is predefined. 

Swarming 

For the bacteria to reach at the richest food location, it 
is desired that the optimum bacterium till a point of 
time in the search period should try to attract other 
bacteria so that together they converge at the desired 
location more rapidly. To achieve this, a penalty 
function based upon the relative distances of each 
bacterium from the fittest bacterium till that search 
duration, is added to the original cost function. Finally, 
when all the bacteria have merged into the solution 
point, this penalty function becomes zero. The effect of 
swarming is to make the bacteria congregate into 
groups and move as concentric patterns with high 
bacterial density. 



(lg) Reproduction 



The original set of bacteria, after getting evolved 
through several chemo tactic stages reaches the 
reproduction stage. Here, best set of bacteria gets 
divided into two groups. The healthier half replaces 
with the other half of bacteria, which gets eliminated, 
owing to their poorer foraging abilities. This makes the 
population of bacteria constant in the evolution 
process. 

Elimination and Dispersal 

In the evolution process, a sudden unforeseen event 
can occur, which may drastically alter the smooth 
process of evolution and cause the elimination of the 
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set of bacteria and/or disperse them to a new 
environment. Most ironically, instead of disturbing the 
usual chemotactic growth of the set of bacteria, this 
unknown event may place a newer set of bacteria 
nearer to the food location. From a broad perspective, 
elimination, and dispersal are parts of the population 
level long distance motile behavior. In its application 
to optimization, it helps in reducing the behavior of 
stagnation often seen in such parallel search algorithms. 
The detailed mathematical derivations as well as 
theoretical aspect of this new concept are presented in 
(Fogel, 1995; Ali and Abd-Elazim, 2011). Fig. 4.shows 
the flow chart of BFOA algorithm. 

Bacterial Foraging Algorithm 

The algorithm of this technique involves two steps. 
[Step 1] Initialization 

i) p is the number of parameters to be optimized. 

ii) S is the number of bacteria to be used for searching 
the total region. 

iii) Ng is the swimming length after which tumbling 
of bacteria will be undertaken in a chemotactic loop. 

iv) is the number of iteration to be undertaken in a 
chemotactic loop ( Nq > Ng ). 

v) N re is the maximum number of reproduction to be 
undertaken. 

vi) N e j is the maximum number of elimination and 
dispersal events to be imposed over the bacteria. 

vii) P ^ is the probability with which the elimination 
and dispersal will continue. 

viii) P (1-p, 1-S, 1) is the location of each bacterium 
which is specified by random numbers on [-1, 1]. 

ix) The value of C (;') which is assumed to be constant 
in this case for all the bacteria to simplify the design 
strategy. 

x) The values of d attract , co amact , h repdmt and a> repelent ■ 

[Step-2] Iterative algorithm for optimization 

This section models the bacterial population 
chemotaxis, swarming, reproduction, elimination and 
dispersal (initially, j=k=l=0). For the algorithm 

updating l automatically results in updating of P. 
[1] Elimination-dispersal loop: 1=1+1 
[2] Reproduction loop: k=k+l 



[3] Chemotaxis loop: j=j+l 

a) For i=l, 2... S, calculate cost function value for each 
bacterium ias follows. 

■ Compute value of cost function J(i, j,k,l) . 

Let J s JiJ,kJ) = J(i,j,kJ) + J cc (O i (j,kJ),P(j,kJ)). 
J cc Is defined by the following equation 

S i 
J cc (0,P(j,k,D))= X J cc (0,0 l (j,k,i)) 

i = 1 



S 

= I 
i = l 

S 

+ I 
i=\ 



-d 



attract 



exp 



P 

Z( 

m = l 



- a attract ^m' ^ 2 



Kepelent exp 



1 ; ? 

-co (o —0 i 

repelent ^ ^ m m' 
m = l 



(21) 



^ et ^last = ^sw^' ^' ^ to save ^ s va ^ ue smce 
one may find a better cost via a run. 

End of For loop 

b) For i=l, 2,....S take the tumbling/swimming 
decision 

Tumble: generate a random vector 

A(;') e S R P with each element A m (i) m=l,2,...p, 



Move: Let 



i U + l,k,D = i U,k,D + C(i) 



Hi) 



a/a 7 (0A(0 



Fixed step size in the direction of tumble for bacterium 
i is considered. 

Compute J(i, j + 1, k, I) and 

J m {iJ + lXl) = J{iJ + lXl) + J cc {0\j + l,k,l),P{j + l,k,l)) 
Swim 

i) Let m=0 (counter for swim length). 

ii) While m<N $ (have not climbed down too long) 

• Let m=m+l 

• If J sw (i,j + l,k,l)<J [ast (if doing better), let 
J te ,=' w (' - .J+U,Z)andlet 

A(/) 



l (j + l,k,l) = l (j,k,l) + C(i) 



^A T (i)A(i) 



and use this l ( j + l,k,l) to compute the 

new J(i,j + l,k, I) 
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• Else, let m = N s . This is the end of the while 
statement. 

iii) Go to next bacterium (i+1) if i it S 

[4] If )<N C , go to [step 3]. In this case, continue 
chemotaxis, since the life of the bacteria is not over. 

[5] Reproduction 

a) For the given k and 1, and for each i=l,2, ..S, let 



•' 7, 



health 



mm \J (i, i,k,l)\. ,, , ,„ 
it sw K,J ' J >be the health 

1..JV 



b) 



of the bacterium i( a measure of how many 
nutrients it got over its life time and how 
successful it was at avoiding noxious 
substance). Sort bacteria in order of ascending 

cost health ■ 

S r =57 2 bacteria with highest 

health va ^ ues ^ e other S bacteria with 



The 
./ 

the best value split. 

[6] If k<N re , go to [step 2]. In this case, one has not 

reached the number of specified reproduction steps, so 
one starts the next generation in the chemotactic loop. 

[7] Elimination-dispersal: for i = 1,2,...., N, with 
probability P ^ , eliminate and disperse each 

bacterium, and this result in keeping the number of 
bacteria in the population constant. To do these, if you 
eliminate a bacterium, simply disperse one to a 
random location on the optimization domain. If 
1< N e j , then go to [step 2]; otherwise end. 

The detailed mathematical derivations as well as 
theoretical aspect of this new concept are presented in 
(Passino, 2002; Mishra , 2005; Fogel, 1995). 

Overview ofPSO 

PSO is a stochastic optimization technique that draws 
inspiration from the behavior of a flock of birds or the 
collective intelligence of a group of social insects with 
limited individual capabilities. In PSO a population of 

particles is initialized with random positions X. and 
— > 

velocities y. , and a fitness function using the particle's 

positional coordinates as input values. Positions and 
velocities are adjusted, and the function is evaluated 
with the new coordinates at each time step. The 
velocity and position update equations for the d-th 
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Initialize all variables. Set all loop counters 
and bacterium index i equal to 0. 



Increase elimination and dispersal loop 
counter 1 = 1+1 



Print the results 
and stop 




Increase reproduction loop counter k = k + 1 




Perform elimination 
and dispersal for all 
bacteria 



Increase chemotaxis loop counter j=j + 1 ^ 



Perform 
reproduction 




^ Increase bacterium index i = i + 1 




Yes 



Compute the cost function value for ith bacterium as 
J ( i, j, k, l ) adding the cell to cell attractant effect to 



nutrient concentration and set J_last = J 



Tumble (let the ith bacterium take the step of height 
C(i) along a randomly generated tumble vector A(i)) 



Compute the objective function value J(y + i,k,i) 
taking in to account the cell to cell attractant effect 



Set swim count m=0 




Set J_last = J(i,j+i,k,i) and swim (let the ith bacterium take a set 
of height C(i) along the direction of the same tumble vector 



FIG. 4 FLOW CHART OF BFOA 
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dimension of the i-th particle in the swarm may be 
given as follows: 



(22) 
(23) 



V' +1 > = *V'> + W J /y- J H<'» + 
x id« + v = x id^ +v id( t + v 

Where 

^id^ t + ^' ^id(t) '■ The updated and current particles 
velocities, respectively. 

X id( t + ^' X id(t) : The updated and current particles 
positions, respectively. 

X lid ' x gd : "^ ae l° ca l anc ^ global best position of each 
particle. 

Cj,C 2 : are two positive constants. 

<Py,<P2 '■ are unit random numbers within the range [0, 1] 

w: is the inertia weight. 

Fig. 5. shows the flow chart of PSO algorithm. 



Start 



Set number of particle and maximum iteration 



Initial population with random position (x) 
and velocity (v) 



Evaluate fitness function for each particle 



Next iteration 



Search for local best position of each 
particle& Search for the global 



Update position 




Update velocity 
J 



Gbest= parameters of best solution | 



0"0 



Fig. 5 FLOW CHART OF PSO ALGORITHM 



Overview ofBSO 

BSO combines both algorithms BFOA and PSO thus 
using advantages of both techniques. The aim is to 
make use of PSO ability to exchange social information 
and BFOA ability in finding a new solution by 
elimination and dispersal. In BFOA, a unit length 
direction of tumble behavior is randomly generated 



which may lead to delay in reaching the global 
solution. In the BSO, technique the unit length random 
direction of tumble behavior can be obtained by the 
global best position and the best position of each 
bacterium by PSO algorithm. The proposed BSO 
algorithm to search optimal values of parameters is 
shown in Fig. 6 and described as follows (Korani, 2008; 
Biswas, Dasgupta, Das, and Abraham, 2007;). 

[Step 1] Initialize parameters 

n, S, N c ,N re , N ed ,P ed , C(i)(i = 12, ,N)J. 

[Step 2] Update the following 

J (i, j,k) : Cost or fitness value of the i-th bacterium 
in the j-thchemotaxis, and the k-th reproduction 
loop. 

-> 

6 g _best : Position vector of the best position 
found by all bacteria. 

J bes( (i,j,k): Fitness value of the best position 

found so far. 
[Step 3] Reproduction loop: k = k + 1 
[Step 4] Chemotaxis loop: j = j + 1 

[Sub step a] For i=l, 2,..., S, take a chemotaxis step 
for bacterium i as follows. 

[Sub step b] Compute fitness function, J(i,j,k). 

• Let m = m + 1 

• If J(i, j + 1, k) < J lasf (if doing better), 

• Let = J(i, j + l,k) and let 



0(i,j + l,k) = 0(i,j,k) + C(i) 



A(0 



a/a t 0-)A(0 



and 



use this 6{i, j + 1, k) to compute the new 
J(i, j + 1, k) as shown in new [sub step f] 

• Else, let m = N s . This is the end of the while 
statement. 

[Step 5] Mutation with PSO operator 

For i=l, 2, , S 

-> 

• Update the 6 g _best and J^ gst 

• Update the position and velocity of the d-th 

coordinate of the i-th bacterium according to 
the following rule: 
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Start 




Initialize all variables. Set all loop counter and 
bacterium index i to 

J 



Increase elimination and dispersal loop 
counter 1 = 1+1 



Print the 
results and 
stop 




Increase reproduction loop counter k = k + 1 




Perform elimination 
and dispersal for all 
bacteria 



Increase chemotaxis loop counter j=j + 1 



Perform 
reproduction 




^ Increase bacterium index i = i + 1 




Yes 



Compute the cost functionvalue for ith bacterium as J ( L 

j, k , i ) and set JJast = J ( ij, k , i j 



Update the direction and position of i-th 
bacterium by PSO 



1 Tamble j 



Compute the cost function value J ( i; k _ i 



'I Set swim counter m = 

Yes 




Swim 



FIG. 6 FLOW CHART OF BSO ALGORITHM 



v id =coV id +c i-n 



6 ~ d 



e n d ew {i,j + \,k) = e° d ld {i,j + \,k) + v™ w 



[Step 6] Let S _r = SI2 

The S _ r bacteria with highest cost function (/) values 
die and other half bacteria population with the best 
values split. 

[Step 7] If k < N re , go to [step 1]. One has not reached 
the number of specified reproduction steps, so one 
starts the next generation in the chemotaxis loop. 

More details of BFOA and PSO parameters are 
presented in Appendix. 

Results and Simulations 

In this section, the superiority of the proposed BSO 
algorithm in designing SSSC (BSOSSSC) in compare to 
optimized SSSC with PSO (PSOSSSC) and optimized 
SSSC controller based on BFOA (BFSSSC) is illustrated. 

Fig. 7.shows the variations of objective function with 
different optimization techniques. The algorithm is run 
keeping limiting value of cost function at 10 5 . The 
objective functions decrease monotonically over 
generations of BFOA, PSO and BSO. Moreover, BSO 
converges at a faster rate (43 generations) compared to 
that for PSO (62 generations) and BFOA (75 
generations). 

Computational time (CPU) of all algorithms is 
compared based on the average CPU time taken to 
converge the solution. The average CPU for BSO is 
21.34 second while it is 27.42 and 34.26 second for PSO 
and BFOA respectively. It is clear that, average 
convergence time for BSO is less than other methods. 
The higher computational time forBFOA is expected 
due to its dependence on random search directions 
which may lead to delay in reaching the global 
solution. On the other hand, PSO suffers from the 
partial optimism, which causes the less exact at the 
regulation of its speed and the direction. In addition, 
the algorithm pains from slow convergence in refined 
search stage, weak local search ability and algorithm 
may lead to possible entrapment in local minimum 
solutions. 

Table 1, shows the system eigenvalues, and damping 
ratio of mechanical mode with three different loading 
conditions. It is clear that the BSOSSSC shift 
substantially the electromechanical mode eigenvalues 
to the left of the S-plane and the values of the damping 
factors with the proposed BSOSSSC are significantly 
improved to be ( a =-2.291,-1.8016,-1.11) for light, 
normal, and heavy loading respectively. Also, the 
damping ratios corresponding to BSOSSSC controllers 
are almost greater than those corresponding to 
PSOSSSC and BFSSSC ones. Hence, compared to 
BFSSSC and PSOSSSC, BSOSSSC greatly enhances the 
system stability and improves the damping 
characteristics of electromechanical modes. 
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Table 2, shows the controller parameters of DC voltage 
regulator, AC voltage regulator, and the additional 
controller obtained by various algorithms. 
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FIG. 7 VARIATIONS OF OBJECTIVE FUNCTION 

TABLE 1 MECHANICAL MODES AND DAMPING RATIOS FOR 
VARIOUS CONTROLLERS AND OPERATING CONDITION 





BFSSSC 


PSOSSSC 


BSOSSSC 


Light 
load 


-1.4487 ± 7.804j 
(0.183) 


-1.778 ± 8.89j 
(0.196) 


-2.291 ± 10.36) 
(0.216) 


Normal 
load 


-1.222 ± 6.1608j 
(0.195) 


-1.474 ± 6.544j 
(0.22) 


-1.8016 ± 6.915j 
(0.2521) 


Heavy 
load 


-0.656 ± 4.95) 
(0.1314) 


-0.785 + 5.26j 
(0.148) 


-1.11 ± 6.18j 
(0.1768) 



TABLE 2 THE CONTROLLER PARAMETERS FOR VARIOUS 
CONTROLLERS 





BFSSSC 


PSOSSSC 


BSOSSSC 


*Pdc 


0.5398 


0.8594 


1.8077 


Ki dc 


1.7412 


0.1041 


1.814 


K Pac 


1.9673 


2.5581 


3.6947 


Ki 

ac 


0.0527 


1.2907 


0.6113 


K 


0.0994 


0.1432 


0.1016 


T l 


0.761 


0.4984 


0.5655 


T 3 


0.5569 


0.3121 


0.4937 



Response under normal load condition: 

The effectiveness of the performance under 0.2 step 
increase in mechanical torque is applied. Fig. 8 shows 
the response of speed for normal loading condition. 
This figure indicates the capability of the BSOSSSC in 
reducing the settling time and damping power system 
oscillations. Moreover, the mean settling time of these 
oscillations is 3.5, 3.8, and 4.2 second for BSOSSSC, 
PSOSSSC, and BFSSSC respectively. In addition, the 
proposed BSOSSSC outperforms and outlasts 
PSOSSSC and BFSSSC controller in damping 
oscillations effectively and reducing settling time. 




FIG. 8 CHANGE IN SPEED FOR NORMAL LOAD CONDITION 

Response under heavy load condition: 

Fig. 9 shows the system response at heavy loading 
condition with fixing the controller parameters. From 
this figure, it can be seen that the response with the 
proposed BSOSSSC shows good damping 
characteristics to low frequency oscillations and the 
system is more quickly stabilized than PSOSSSC and 
BFSSSC. The mean settling time of oscillation is 3.2, 3.8, 
and 4.1 second for BSOSSSC, PSOSSSC, and BFSSSC 
respectively. Hence, the proposed BSOSSSC extend the 
power system stability limit. 




FIG. 9 CHANGE IN SPEED FOR HEAVY LOAD CONDITION 
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Conclusions 

In this study, a new optimization algorithm known as 
BSO, which synergistically couples the BFOA with the 
PSO for optimal designing of SSSC damping controller 
is thoroughly investigated. For the proposed controller 
design problem, an integral time absolute error of the 
speed, DC voltage and transmission line power of 
SSSC is taken as the objective function to improve the 
system response in terms of the settling time and 
overshoots. The superiority of this objective function is 
that the effect of SSSC signal is taken into 
consideration. Simulation results are presented for 
various loading conditions to verify the effectiveness 
of the proposed controller design approach. Moreover, 
the proposed control scheme is robust, simple to 
implement, yet is valid over a wide range of operating 
conditions. 
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APPENDIX 

The system data are as shown below: 

a) Synchronous generator (p.u) X ^ =1.07, X =1.0, X ^ =0.3, 

T do =5 - 9 ' 11=137 ' p e=0.9, Vt=1.0. 

b) Excitation system K . =400 and T. =0.05 sec. 

c) Transmission line (p.u), X { ^=03, X ^ =0.3 

d) SSSC parameters (p.u), X^ =0.05, \ 'jyQ =1-0, Cj^q =1-0. 

e) Bacteria parameters: Number of bacteria =10; number of 
chemotatic steps =10; number of elimination and dispersal 
events = 2; number of reproduction steps = 4; probability of 
elimination and dispersal = 0.25; the values of ^ a ff raC f =0.01; 

the values of 03^^ =0.04; the values of =0.01; 
the values of a> repelent =™. 

f) PSO parameters: Cj = =2.0, CO =0.9. 
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